Flow-induced pruning of branched systems 
and brittle reconfiguration 
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Whereas most plants are flexible structures that undergo large deformations under flow, another 
process can occur when the plant is broken by heavy fluid-loading. We investigate here the mecha- 
nism of such possible breakage, focusing on the flow-induced pruning that can be observed in plants 
or aquatic vegetation when parts of the structure break under flow. By computation on an actual 
tree geometry, a 20-yr-old walnut tree (Juglans Regia L.) and comparison with simple models, we 
analyze the influence of geometrical and physical parameters on the occurrence of branch breakage 
and on the successive breaking events occurring in a tree-like structure when the flow velocity is 
increased. We show that both the branching pattern and the slenderness exponent, defining the 
branch taper, play a major role in the breakage scenario. We identify a criterion for branch breakage 
to occur before breakage of the trunk. In that case, we show that the successive breakage of periph- 
eral branches allows the plant to sustain higher flow forces. This mechanism is therefore similar to 
elastic reconfiguration, and can be seen as a second strategy to overcome critical events, possibly a 
widespread solution in plants and benthic organisms. 
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I. INTRODUCTION 

Most living systems are surrounded by a fluid, be it air 
or water. When this fluid flows, it generates mechanical 
forces, that may have major consequences on growth as 
well as on reproduction or survival Typical cases 

are trees subjected to wind or corals subjected to water 
currents. In terms of flow-induced deformations, two typ- 
ical behaviors can be pointed out. In the most common 
one, the solid undergoes large elastic deformations, for in- 
stance in crops or aquatic vegetation. In the second type, 
the system breaks before any significant deformation can 
occur; this will be referred to as brittle behavior in the 
following. The former has been abundantly studied, a 
key result being that of load reduction by elastic recon- 
figuration [4, rj. The latter has already been described in 
trees or corals J(| 0] , but to the best of our knowledge the 
effect of branching has never been studied theoretically. 
Therefore, we shall focus hereafter on brittle branched 
slender systems, which are ubiquitous in nature: trees 
H, bushes, algae @, corals @ and corallines [ioj . to list 
a few. In the following we refer mainly to trees under 
wind loading, with the understanding that these results 
are also applicable to a large variety of other biological 
systems under fluid-loading. 

For a brittle branched system attached to a support, 
breakage under flow may occur in three distinct types: 
(i) base breakage, Fig. fTJi, when the attachment to the 
ground is broken, as in uprooting, (ii) trunk breakage, 
Fig-Hb, when the main element is broken, and (iii) branch 
breakage, Fig. [Tp, when an upper element breaks, as in 
flow-induced pruning. 
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In fact, the distinction between trunk and branch 
breakage has a biological relevance, since breakage of the 
trunk is likely to be fatal, while re-growth is often pos- 
sible after branch breakage. Moreover branch breakage 
does reduce loads on the trunk and the attachment, as 
in elastic reconfiguration, thereby delaying their break- 
age @, [ll[ . Finally, branch breakage can also be part of 
the asexual reproduction process by propagation. This is 
observed in terrestrial plants such as willows and poplars 
[l^ . and in stony corals such as Acropora Cervicornis or 
Acropora Palmata [13l. Il4j. 

Breakage is the consequence of an unacceptable stress 
level; it is therefore directly related to the stress state in 
the structure [ll|, [HI- In particular, the issue of whether 
the stress level is uniform or not in the tree is crucial, 
as breakage is expected to occur at the point of maxi- 
mal stress. For instance, Niklas and Spatz [ll| showed 
that in a cherry tree the stress level varies by one to 
two orders of magnitude within the tree and has a local 
maximum in the branches. On the other hand, Bejan 
et al. showed that the flow-induced stress is uniform for 
a tapered trunk when the taper is linear (l6| . In fact 
the stem taper is an important parameter regarding the 
stress distribution; see the discussion in 

Several questions remain however regarding the flow- 
induced breakage of tree-like structures: (i) what are the 
effects of the geometrical and physical parameters on the 
occurrence of branch breakage? (ii) How do the breaking 
events occur successively as the flow is increased? (iii) 
Assuming that branch breakage is favorable in biologi- 
cal terms, is it compatible with other constraints on the 
geometry? The aim of this paper is to address these ques- 
tions, using simple numerical and analytical models for 
the mechanical behavior of slender and brittle structures. 
The modeling assumptions and framework used through- 
out the paper are first presented in Section [TXJ In Section 
IIII1 we compute the stress distribution and successive 
breaking events in a complex tree, using the geometry of 
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an actual walnut tree. Using an idealized branched sys- 
tem, we derive conditions for branch breakage in Section 
IIV1 These are further analyzed for a tapered beam, here 
referred to as the slender cone model, in Section [V] The 
corresponding three geometries are sketched in Fig. [2] 
Finally a general discussion and conclusion are given in 
Section Ell 

II. MECHANICAL MODEL AND 
PARAMETERS 

Throughout the paper, we consider a cross flow over 
the entire structure, uniform, as the dependence of the 
stress on the wind velocity profile was shown to be small 
[TT| . Also, only static loads are taken into account, and 
the corresponding fluid force magnitude / per unit length 
reads 

/ = \ P C D DU\ (1) 

where U is the free stream velocity, p its density, D the 
local branch diameter and Cd the drag coefficient @, [2(| ■ 
The direction is assumed to be that of the flow velocity. 
The fluid load is here computed on a leafless branch, and 
the influence of leaves will be discussed in Section IVI1 

This load is applied on the whole branched system, 
which is held by a perfect clamping at the base. Because 
of the high slenderness of the system, we use a standard 
linear beam theory to derive the stress state, essentially 
the bending moment M. The maximum stress in the 
cross-section resulting from this bending moment is the 
skin stress, defined as £ = 32M/ttD 3 [2lll23|. 

The brittle behavior is introduced as follows: (i) the 
deformations are assumed to be negligible, so the stress 
state is computed on the initial configuration, without 
elastic reconfiguration, (ii) when increasing the flow ve- 
locity U, breakage occurs when and where the local skin 
stress £ reaches a critical value, £ c . Then, the broken 
branch is removed, and this results in a new flow-induced 
stress state. Flow velocity may then be further increased 
until a new breaking event occurs. 



Throughout the paper, the relevant dimensionless 
number to scale the fluid-loading pCnU 2 with respect 
to the critical stress £ c is the Cauchy number, defined as 

Cy = &^G, (2) 

where G is a geometrical factor introduced for compar- 
ison purpose and defined such that £ — £ c at the base 
of the intact structure when Cy = 1- Note that this 
Cauchy number is similar in principle but differs from 
that used in the analysis of flow-induced elastic deforma- 
tion, namely Cy = pC D U 2 /E @, [|; the critical stress 
£ c simply replaces here the Young modulus E. 

The non-dimensional stress is defined as a — £/£ c and 
the non-dimensional bending moment as m = M/M c , 
with M c — Ti c nDg/32, Db being the base diameter [22| . 
This latter scaling is chosen so that failure occurs at the 
base of the trunk when m = 1. The non-dimensional 
vertical coordinate z is defined using H, the height of 
the structure, as a reference length scale. 

III. FLOW-INDUCED PRUNING OF A 
WALNUT TREE 

The geometry of the branched system is expected to 
have a large influence on the stress state and thus on the 
location and timing of breaking events. We therefore first 
apply the procedure described above using the digitized 
geometry of an actual 20-yr-old walnut tree (Juglans Re- 
gia L.) described in (Fig. Oi). This tree is 7.9 m 
high, 18 cm in diameter at breast height (dbh), and has a 
sympodial branching pattern [2 31 ] and about eight orders 
of branching. The stress state under flow is computed 
using a standard finite element software (CASTEM v. 
3M [24|). and is presented in Fig. ^jp for four different 
branching paths. 

We observe that the stress level is not uniform but 
shows a maximum located in the branches, which is con- 
sistent with the results of Niklas and Spatz |Tl| which 
are sketched in Fig. [3^,. Note that since a varies linearly 
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FIG. 2: Geometries of the models used in the paper: fa) Section IIIII Walnut tree, as in [18J; (b) Section HVl Idealized branched 
system, as in (c) Section [V] Tapered beam, as in [a.[lr3|. 




(a) cr (b) cr (c) 

FIG. 3: Non-dimensional stress profile cr in a tree under cross-flow, (a) Schematic view of the stress profiles given by Niklas and 
Spatz f f] for cherry trees, showing a local maximum near the top. (b) Computed stress profiles along four branching paths, 
A (x), B (□), C (A) and D (o) in the digitized tree geometry shown in (c). 



with the fluid-loading Cy, one needs only to focus on 
the critical situation where cr = 1 is first reached in the 
structure. In this tree, the criterion for breakage is sat- 
isfied first in a branch and not in the trunk. This corre- 
sponds to the mechanism of branch breakage, as denned 
in Section HI If the fluid- loading is further increased after 
removal of the broken parts, successive breaking events 
are observed, in a flow-induced pruning sequence: Fig.0^ 
shows three states of the tree at increasing Cauchy num- 
ber with branches progressively removed as they break 
off. 

During the sequence of breakage, the bending moment 
at the base of the tree, mb, evolves significantly with the 
Cauchy number, Fig. @]d. Up to the first breakage, the 
moment is proportional to the fluid-loading Cy (zone I 
in Fig. \%jp) . Then, in a small range of load increase (zone 



II) , all large branches are broken at an intermediate level, 
resulting in a significant decrease of the bending moment. 
Breakage then continues but to a much smaller extent 
(zone III) , while the moment increases almost linearly up 
to the value m-b = 1 when the trunk breaks. Note that 
the benefit of this sequence of breaking events is that the 
critical value of the base moment — 1 is reached only 
at Cy — 10 instead of Cy = 1 if there was no branch 
breakage. This corresponds to more than a factor of 3 on 
the acceptable fluid velocity. For instance, for a critical 
stress S c = 40 MPa, which is the order of magnitude of 
maximum acceptable bending stresses measured in trees 
[l2l [25j , the maximum sustainable fluid velocity before 
trunk breakage is increased from C/ ~ 30 m.s -1 without 
branch breakage to U ~ 100 m.s -1 with branch breakage. 

To summarize, this set of computations clearly shows 
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FIG. 4: Computed sequence of branch breakage in the wal- 
nut tree: (a) A: initial tree for CV < 0.67; B: after breakage 
in large branches, CV = 1.7; C: just before trunk breakage, 
CV = 10.7. (b) Corresponding evolution of the bending mo- 
ment at the base of the tree mb, in three distinct ranges. 
The dashed line shows the moment that would exist without 
breakage. The dotted line shows the critical value m,b that 
causes trunk breakage. 



that branch breakage can occur prior to trunk breakage, 
and that the sequence of flow- induced pruning results in 
a significant reduction in the load applied on the base of 
the tree, or equivalently, an increase in the sustainable 
fluid velocity. To further analyze this process, we turn to 
a simple model in the next section. 



IV. THE IDEAL TREE MODEL 

A. Infinite branched tree 

To establish the relation between the parameters of the 
system and the flow-induced pruning process, we simplify 
the problem to its essential elements: the branched geom- 
etry and the slenderness of branches; we disregard here 
the effect of branch orientation relative to the flow. Sim- 
ilarly to [H, we consider first an infinitely iterated sym- 
podial tree made of cylindrical branches (Fig. [5]). Two 
parameters only are needed to describe this ideal tree: (i) 
the branching ratio A, giving the reduction of diameter 
through branching, and (ii) the slenderness exponent /3, 
giving the relationship between length and diameter in 



FIG. 5: Idealized branched system, (a) Infinite iterated tree. 
The sub-tree II is equivalent to the whole tree I but for a 
change of scales, (b) Finite iterated tree and corresponding 
notations. 



branch segments of the tree, so that 
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where and Lfc are the diameter and length of a branch 
segment of order k, see Fig. [S^, [19j. Typical values of 
these parameters are A < 1 and 1 < (3 < 2. Note that 
the number of branches emerging from a branching point 
is typically equal to 1 / A (2(| . 

We use now a scaling argument similar to that of 
19] for the dynamics of trees. On the ideal infinitely 
branched system of Fig. [5ji, we can compare the stress 
level in branch k = 1 (the trunk) and in branch k = 2. 
The sub-tree labeled II in Fig. [S^ is identical to the full 
tree, I, but for a change in length and diameter scales. 
All diameters (resp. lengths) in II are reduced by a fac- 
tor A 1 / 2 (resp. A 1 / 2 ' 3 ). Let Si be the maximum skin 
stress in the trunk (fc = 1) under a given fluid-loading 
U, and £2 the maximum skin stress in the branch k = 2. 
The relations between the flow velocity and £1 or £2 
are identical, but for the change of diameter and length 
scales. The dependence of the stress on diameter and 
length is the following: (i) £ varies as M/D 3 , where M 
is the bending moment, (ii) M varies as fL 2 , where / is 
the norm of the local fluid force, Eq. ((T]), (iii) / varies as 
pU 2 D. Hence £ varies as pU 2 (L/D) 2 . We therefore may 
state that 



£2 
£1 



D 2 



£1 



= A" 



(4) 



Since A < 1, the condition for the stress to be higher in 
branches than in the trunk becomes 



/? > 1. 



(5) 



Here the only parameter controlling the possibility of 
branch breakage is the slenderness exponent, a classical 
parameter in the allomctry of trees. As /3 is typically 
greater than 1 for trees, branch breakage is expected 
to occur. This simplistic approach now deserves to be 
improved, as the assumption of an infinite number of 
branching levels is very strong, and may not be compati- 
ble with the constraint that the tree area has to be finite. 
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B. Finite branched tree 

Let us consider now the same idealized tree, but with 
a finite number of branching iterations (Fig. [SJd). This 
structure has ./V levels, which are labeled in this section 
from the top to the bottom. Note that n — N — k + 1 , 
where n is the label of the previous section from the base 
of the tree. The trunk corresponds now to the last level, 
N. At each level n, we define the branch diameter D n 
and length L n , which can be expressed as a function of 
the trunk diameter and length Dn and Ln as 



N—.. 



N, 



L, 



N- 

A 2 <3 L 



N- 



(6) 



By a simple integration of the fluid force on the 
branches, the moment at the base of a branch of order n 
may be derived, as well as the corresponding skin stress, 
which is obtained in non-dimensional form as 



B\% +CAV n 



CyA^ JV ( 
where the Cauchy number Cy is defined as 

pC D U 2 



C Y = 



Ln 
D N 



(7) 



(8) 



and A, B and C are functions of (3 and A only. The 
detailed derivation of Eq. (0 as well as the expression of 
A, B and C can be found inlA"! 

A systematic numerical exploration of the (A, /3) pa- 
rameter space shows that when j3 < 1 the stress always 
increases from top to bottom. Conversely, for j3 > 1 , the 
stress reaches a maximum at branch level n c and then 
decreases from top to bottom, provided that N > n c , 
where n c depends on A and (3. This dependence is given 
in Fig. [6] This analysis with a finite tree model gives a 
criterion consistent with that of the infinite tree model, 
namely (3 > 1. Moreover, the other parameter, A, is 
found to affect only the location of possible breakage. 
This suggests that branching is not a key factor in the 
occurrence of branch or trunk breakage. In the next sec- 
tion we explore a simpler model of the slenderness effect. 



V. THE SLENDER CONE MODEL 
A. Flow-induced stress 

The simplest model that allows one to take into ac- 
count a relation between diameters and lengths through 
a slenderness exponent is a cone. This formulation is re- 
lated to MacMahon and Kronauer's equivalent geometry 
of a tree, a tapered beam with a rectangular cross-section 
of dimensions varying as power laws of height [1, [13] ■ 

The geometry considered here is a slender cone with 
a circular cross-section, Fig. [7^,, and we follow the same 
mechanical approach as for the previous geometries. Let 
H be the cone height, dn = L>h/H the dimensionless 
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FIG. 6: Location of the maximum of stress under cross-flow 
in an idealized tree model, as a function of the slenderness 
exponent (3 and the branching parameter A. The location is 
given in the form of the number of branching levels counted 
from the top of the tree, Fig. [5b. For f3 < 1, the breakage is 
directly at the base of trunk. 



diameter at the base and z the vertical coordinate which 
is orientated downwards in this section. The cone dimen- 
sionless diameter is given by 



d(z) = d H z?. 



(9) 



Using the same formulation as in the previous section, 
the stress state along the cone is obtained as 



where the Cauchy number is defined here as 
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pC D U 2 



(10) 



(11) 



From Eq. (fT0|) . we readily observe that: (i) for (3 = 1, 
the constant stress case of Bejan et al. [l6| is found; (ii) 
for j3 < 1 the stress increases with z and is therefore 
maximum at the base, Fig. [7^; (iii) for j3 > 1 the stress 
decreases with z, and the maximum, discussed further, 
is not at the base, Fig.[7J>c. These results are consistent 
with the condition for branch breakage in the previous 
section. 

To avoid the singular case of infinite stress at z — for 
(3 > 1, we use a cone truncated at z — zo, Fig. [Th. The 
truncation zq corresponds to the first breakage occurring 
as soon as U ^ 0, and its value is chosen arbitrarily. The 
corresponding stress state is then 



o(z) 
C Y 



z 2(i-/3)_ (2 + / ^i 



i+V- 



3/J 



+ (l + (3)z 2 ^z'^, 



(12) 

which reduces to Eq. (fTUf when zq = 0. The detailed 
derivation of this equation is given in|B] For (3 > 1, the 
stress shows a maximum before decreasing downwards, 
as illustrated in Fig. [7fc. The limit case Zo = is in 
fact equivalent, in the ideal tree model of Section HVl to 
the limit as N goes towards infinity, which would lead 
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(c) 

FIG. 7: The slender cone model: geometry and stress profile 
under uniform cross flow, (a) cone with f3 < 1 (here 0.75), 
showing a maximum of stress at the base; (b) cone with f3 > 1 
(here 2), showing a maximum at the top; (c) cone truncated 
arbitrarily at zq — 0.3 showing a local maximum. 



to a vanishing diameter at the tip. There is therefore an 
analogy between the cone truncation and the ideal tree 
with a finite number of branching levels. 



B. Sequence of breaking events 

Considering now the generic case of the truncated cone, 
Fig. [7J;, we analyze the sequence of breaking events re- 
sulting from an increasing fluid- loading Cy ■ The stress a 
increases linearly with Cy up to the point where its max- 
imum value reaches the limit of breakage, a = 1. This 
defines the first breaking event at Cy = C Y occurring 
at z = Z\. It results in a new truncated cone, and the 
process is repeated as Cy is further increased. Eventu- 
ally, when the cone becomes truncated close to the base, 
the maximum stress may be reached at the base itself, 
resulting finally in base breakage. 

This sequence of breaking events may be analyzed in 
terms of the maximum fluid-loading Cy ax that the cone 
can support before breaking at the base. As illustrated in 
Fig. [51 this is strongly dependent on (3. When (3 < 1, the 
first breaking event is at the base so that C™ ax = 1 . Con- 
versely when (3 > 1, breaking occurs progressively as Cy 
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FIG. 8: Maximum fluid load that the cone can support as a 
function of the slenderness exponent. Note that for ft > 1 the 
curve is the lower bound of all possible evolutions. 



is increased, and the base breakage is delayed, C™ ax > 1. 
The precise value of Cy where the base breaks depends 
on the initial truncation zq, but is always higher than a 
lower bound that can be computed from Eq. (| 1 2[) . which 
is shown in Fig. |8j We observe a significant increase of 
the ability of the system to sustain fluid-loading when 
/3>1. 

In terms of base moment, the sequence of breaking 
events can be easily computed, Fig.El For (3 < 1 the base 
moment increases linearly with Cy until base breakage 
occurs, mb = 1 for Cy = 1. For f3 > 1 the sequence 
of breaking events results in sudden drops in base mo- 
ment followed by linear increase up to the next breaking, 
as illustrated in Fig. [S] Since the sequence of breaking 
events is a discrete process that depends on the initial 
truncation zq, there exists, for a given Cauchy number 
Cy , a wide range of acceptable cone heights and thereby 
a wide range of corresponding base moments. In prac- 
tice, for all possible values of zq, the evolution of mb 
remains bounded between its values for the shortest and 
highest cone that can exist at each Cauchy number. This 
is represented by the shaded region in Fig. [9] 

These results show that the simple cone model contains 
the key elements to understand the effect of geometry on 
(i) the stress profile, (ii) the sequence of breaking events 
and (iii) the consequences on the evolution of base load 
when the fluid velocity is increased. Here again, the es- 
sential criterion concerns the slenderness exponent j3. 



VI. DISCUSSION AND CONCLUSIONS 

Starting from the case of a full walnut tree geometry, 
we have used models of increasing simplicity. This al- 
lowed us to point out the role of various parameters on 
the process of breakage under fluid-loading. The first is- 
sue that had to be addressed was that of the flow-induced 
stress distribution. As noted by other authors, the stress 
is not necessarily maximum at the base pj] . [l6| . In fact in 
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FIG. 9: Moment at the base of the cone as the fluid-loading is 

increased. ( ) direct base breakage occurring when /3 < 1; 

( — ) progressive breaking for /3 > 1 (here ft = 2). The shaded 
region shows all possible values depending on the initial trun- 
cation zq. The cone state is shown for three values of CV- 



the walnut tree of Section UTTl the stress has a local max- 
imum at about mid height. Using the ideal tree model 
in Section IIV1 we have shown that the existence of this 
maximum is related to the value of the slenderness ex- 
ponent, /3, being larger than one: in fact this allometry 
parameter is about 1.37 for this particular walnut tree 
19]. Following Bejan et al. [lfjj], we recover the critical 
value of j3 = 1 in the simplest model, that of a cone in 
Section El 

Actually, some refinement is needed here to understand 
the precise location of the maximum of stress. We have 
shown in Section IIVI that the location of this maximum 
was also dependent on the branching parameter A, in 
the form of the parameter n c , which is the number of 
branching levels from the top to this maximum point. 
For our walnut tree, where A ~ 0.25, we obtain n c = 6 
using Fig. [6] This is smaller than the total number of 
branching levels in the walnut tree which is about 8 [l8[ . 
A local maximum of stress is therefore expected in the 
branches, and is actually observed in Fig. [3J 

The second issue was that of the sequence of breaking 
events occurring when the fluid- loading Cy is increased. 
Using a brittle fracture model for the walnut tree in Sec- 
tion IIII1 we have shown that most large branches broke 
in a short range of flow velocity, and that breakage of 
the trunk occurred much later. The large size of bro- 
ken branches can be explained by the value of n c — 6 
found above. All large branches do not break exactly at 
the same value of the Cauchy number. This is due among 
other reasons to some variability in the allometry param- 
eters A and /3 within the tree. Once all large branches are 
broken, the remaining tree shape, C in Fig. does not 
have enough branching levels to have a local maximum, 
and the next breaking event occurs at the base of the 
trunk. Note that the process of branch breakage in the 
walnut tree allowed the tree to have a much larger accept- 
able Cauchy number before breakage of the trunk. This 



can also be analyzed using the cone model as in Section 
fVl where the critical Cauchy number for base breakage 
is clearly dependent on (3 (Fig. [5]) . 

The third issue was that of the evolution of the load 
at the base of the tree. For the walnut tree, Fig. @J), 
the sequence of successive breakage of the large branches 
results in a significant decrease of the drag-induced mo- 
ment at the base. This can be understood using the cone 
model, where the sequence of breaking event and corre- 
sponding drops of base moment can be tracked, Fig. |H1 
We may therefore state that the essential characteristics 
of branch breakage and corresponding load evolution in 
the walnut tree can be understood using our simple ideal 
tree model and cone model. 

The analytical results of Sections IIVI and |V] were ob- 
tained considering that all parameters have self-similar 
variations. However, this was not the case for the wal- 
nut tree computations of Section iLUl which suggests that 
the behaviors pointed out in this study can be general- 
ized to structures that do not necessarily have self-similar 
variations of their parameters. Moreover, the ideal tree 
and cone models can be easily extended to incorporate 
other features of the problem, such as a dependence of 
all parameters with z: the flow velocity U, the material 
properties through the critical parameter S c , and even 
the drag coefficient Co, which allows one to take easily 
into account the additional drag of leaves. Preliminary 
results, not shown here for the sake of brevity, showed 
that the criterion for branch breakage takes the same 
form, but involves both /3 and the corresponding param- 
eter related to the additional z-dependence. Taking into 
account a significant elastic deformation before load frac- 
ture, or incorporating dynamical effects, would be much 
more complex. 

Considering the simplicity of the criterion that we have 
found for branch breakage, we can test whether it is gen- 
erally satisfied. MacMahon and Kronauer 27J have noted 
that (5 is usually larger than 1 and typically around 1.5, 
while A is typically close to 0.25. This leads to a max- 
imum stress located at a branching level n c = 5 count- 
ing from top down. This is clearly in the branches as 
trees generally have more than 5 orders of branching. 
We may therefore state that branch breakage can be ex- 
pected in most sympodial trees. This is illustrated in 
Table |TJ where the values of parameters are given for sev- 
eral trees. 

Clearly the possibility of branch breakage is favorable 
in terms of survival of an individual tree in the face of ex- 
treme fluid-loading. It may also be favorable in terms of 
tree development by removing the less vigorous branches. 
The question then arises as to whether this implies new 
constraints on the geometry of the tree. It appears from 
our results that the constraint (3 > 1 is not incompatible 
with other constraints such as the optimal resistance to 
buckling under gravity, which requires j3 = 3/2 Q. The 
same result was obtained considering the wind effect on 
trees but for an overcrowded tree canopy [l7| . Similarly 
/3 > 1 is compatible with a constraint for optimal dis- 
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>nicss Branching Total orders of Predicted branch Predicted 
ent parameter A branching N breakage level n c breakage type 



[18 19] Walnut Tree 



Red Oak 
d± Qaeraus Rabr, 



White Oak 1 
Quercus Alba 



White Oak 2 
Quercus Alba 



Poplar Tree 1.5 
Populus Tremoloides (estimated) 



Pin Cherry 
Linus Pensylvfm 



White Pine 
Pin-as Strobu 



TABLE I: Predicted breakage type using the results of Section HVl Branch breakage is predicted when n c < N. 



sipation [lj| [28| , that modal frequencies have a ratio of 
less than two, requiring that j3 > 1 for A = 0.25. 

The particular case of branched corals [H, [l3|, EH is 
somewhat different. The segments are similar in length 
and diameter, so that A ~ 1 and j3 ~ 1 in our vari- 
ables, but with a number of branches emerging from one 
branching not equal to 1/A. An analysis similar to that 
of Section ITVl shows that breakage is expected at the bot- 
tom. This is the case in most isolated corals. 

More generally we may place these results in the over- 
all context of reconfiguration, as introduced by Vogel @ . 
This originally referred to the reduction of loading made 
possible by elastic deformation. For a plant, it is a cru- 
cial mechanism to survive heavy fluid-loading. But plant 
tissues are not all very elastic, and plant parts are not all 
very flexible. Our results on the role of branch breakage 
in reducing loading show that, in parallel with elastic re- 
configuration, there exists a mechanism of brittle recon- 
figuration. There are therefore two distinct strategies to 
overcome critical events. The first is evidently reversible 
in the short term by elasticity. The second is also re- 
versible by re-growth, but only in the long term. Thus 
flow-induced pruning is possibly a widespread mechanism 
in plants or benthic organisms that support heavy load- 
ing by the surrounding fluid environment. 
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Appendix A: Stress derivation in finite branched 
tree model 



In order to compute the stress along the finite ideal 
tree, we introduce /„ the fluid force per unit length at 
level n, f n — \pCoU 2 D n: with the same notations as 
Eq. ([T]). At each level n, we consider two force compo- 
nents: (i) the shear force r n in the flow direction and (ii) 
the bending moment M n in the direction normal to the 
flow. Due to the free condition at the top, To — and 
Mq = 0, and for n > 1 



fn L r , 



■pr n -i, 



(Al) 



M„ 



1 



f n L 2 n + p (M n _i + L n r n -i) , (A2) 



where p is the number of branches emerging from one at 
a branching point (p = 1/A). The non-dimensional stress 
er n at level n reads 

° n = v n3 A 3) 

By integration of Eqs. (|Alj) and (|A2j) . the stress at each 
level can be obtained, 



with 



and 



-N 



A\ iL e ln + B\% +CAV' 



8 (Ln_ 
7t \D N 



pC D U 2 



A 



B = 



C = 



A 2 f> + 1 



(A^ - l) ( 
A^ + 1 



A 2 <3 



, 2=1 

A 2 f 



l) (a* 



hrw- - i) (\^ - i) 



(A4) 
(A5) 

(A6) 

(A7) 
(A8) 
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Appendix B: Stress derivation in the slender cone 
model 



with zq > 0. The local non-dimensional skin stress reads 



The stress state for the slender cone model is obtained 
by direct integration of the fluid force defined in Eq. (fl} , 
using Eq. © for the diameter. The shear force and re- 
sulting bending moment read 



r(z) 



f(z')dz', M(z) 



T(z')dz', (Bf) 



a(z) 



32M(z) 
7rS c d(z) 3 



(B2) 



The integration of these equations give Eq. (!ID) and 
Eq. (TT2)) depending on z . 
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